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I. Introduction 

Aircraft noise has become an increasing concern for commercial airlines. Worldwide 
demand for quieter aircraft is increasing, making the prediction of engine noise 
suppression one of the most important fields of research. 

The Low-Pressure Turbine (LPT) can be an important noise source during the approach 
condition for commercial aircraft. The National Aeronautics and Space Administration 
(NASA), Pratt & Whitney (P&W), and Goodrich Aerostructures (Goodrich) conducted a 
joint program to validate a method for predicting turbine noise attenuation. The method 
includes noise-source estimation, acoustic treatment impedance prediction, and in-duct 
noise propagation analysis. 

Two noise propagation prediction codes, Eversman Finite Element Method (FEM) code 
[1] and the CDUCT-LaRC [2] code, were used in this study to compare the predicted and 
the measured turbine noise attenuation from a static engine test. 

In this paper, the test setup, test configurations and test results are detailed in Section II. 
A description of the input parameters, including estimated noise modal content (in terms 
of acoustic potential), and acoustic treatment impedance values are provided in Section 
III. The prediction-to-test correlation study results are illustrated and discussed in 
Section IV and V for the FEM and the CDUCT-LaRC codes, respectively, and a 
summary of the results is presented in Section VI. 

II. Description of Test Configurations and Results 

Test Description 

Engine noise testing was conducted at Pratt & Whitney’s C-ll test stand in West Palm 
Beach, Florida. Figure 1 shows the far-field microphone array, which consisted of 32 
microphones at a 150-ft radius around the engine. The microphones were spaced in 5- 
degree increments over a range of angles between 5 and 160 degrees from the inlet axis 
of the engine. The engine was a PW4000-94, which has a four-stage low-pressure turbine 
with blade counts of 128, 130, 126 and 128 (ffont-to-back). The PW4000-94 family of 
engines has been in production since 1987. 
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Figure 1 : Far field microphone array at P&W C-l 1 Test Stand. 

Figure 2 shows a cross section of the turbine exhaust duct, which includes acoustic liners 
on the nozzle and plug. 



Figure 2: Turbine exhaust duct geometry and liner locations. Dimensions are normalized 

by the exit diameter of the duct. 

To determine the effect of the acoustic liners on turbine noise suppression, the following 
four configurations were tested: 

Configuration 1 : hardwall nozzle and hardwall plug 
Configuration 2: treated nozzle and hardwall plug 
Configuration 3: hardwall nozzle and treated plug 
Configuration 4: treated nozzle and treated plug. 
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Hardwall Spectra 

Figure 3 shows the sound pressure spectra (SPL) measured for Configuration 1 (hardwall 
nozzle, hardwall plug) at angles of 100, 115 and 130 degrees from the inlet axis, with the 
turbine operating at 62% of its design speed. Corrections were applied to the spectra to 
remove the effect of atmospheric attenuation [4, 5]. The spectra are presented versus 
engine order, which is defined as, 



where 


f = frequency (Hz) 

Q = rotational speed (rev/sec) 


To focus attention on turbine noise, the abscissa has been limited to engine orders 
between 100 and 156. The solid black lines represent the measured engine spectra, and 
the dashed blue lines represent the estimated turbine spectra. The estimated turbine 
spectra are based on Pratt & Whitney’s empirical method for separating turbine noise 
from measured engine noise. 

The sound pressure spectra are characterized by a broad haystack and tones at engine 
orders of 126, 128 and 130. The tones are most apparent at the 100-degree angle. 
However, the tones are barely visible at the 1 15-degree angle, and are not observed at the 
130-degree angle. The tones correspond to the blade-passage frequencies of the low- 
pressure turbine. The haystack is due to frequency scattering that occurs when the turbine 
tones pass through the primary and secondary shear layers [6]. The spectral 
characteristics shown in Figure 3 are qualitatively similar for all configurations and 
speeds presented in this paper. 


ANGLE = 100 deg 


ANGLE = 115 deg 


ANGLE = 130 deg 





Figure 3: Sound-pressure spectra, Configuration 1, 62% design speed. 

Liner Effect 

For each configuration, turbine sound-power spectra (PWL) were determined by 
integrating the turbine sound-pressure spectra (SPL) over the polar array of microphones. 
The integration was limited to angles between 85 and 145 degrees from the inlet axis. 
The turbine overall sound-power level (OAPWL) was determined by integrating the 
turbine sound-power spectra over the range of engine orders shown in Figure 3. 
Consequently, the overall sound-power level represents the power associated with the 
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blade-passage frequencies and haystack produced by all stages of the LPT. Liner 
insertion loss was determined by subtracting the treated overall sound-power level (i.e. 
Configurations 2 through 4) from that of the fully hardwall configuration (i.e. 
Configuration 1). 


Figure 4 shows the turbine sound-power spectra for each configuration with the turbine 
operating at 62% of its design speed. Figure 5 shows the turbine overall sound-power 
level versus percent design speed. In Figure 5, the symbols represent data points, and the 
solid lines represent 3 rd -order polynomial curve fits of the data. As expected, 
Configuration 4 (treated nozzle, treated plug) has the largest insertion loss, which is 
between 6 and 7 dB over the range of speeds shown. The insertion loss of Configuration 
2 (treated nozzle, hardwall plug) is between 4 and 5 dB, and the insertion loss of 
Configuration 3 (hardwall nozzle, treated plug) is approximately 1 dB. These results are 
consistent with the expectation that the insertion loss should be proportional to the 
amount of treated area in the turbine exhaust duct. 


— (1) Hardwall nozzle, hardwall plug 
— (2) Treated nozzle, hardwall plug 
(3) Hardwall nozzle, treated plug 
— (4) Treated nozzle, treated plug 



o (1) Hardwall nozzle, hardwall plug 
v {2) Treated nozzle, hardwall plug 
A (3) Hardwall nozzle, treated plug 
° (4) Treated nozzle, treated plug 



Figure 4: Turbine sound power spectra 
at 62% design speed. 


Figure 5: Turbine overall sound 
power level. 


Table 1 shows the insertion loss at the engine speeds used for the FEM and CDUCT 
analyses. The tabulated insertion loss values are based on the curve fits shown in Figure 
5. 


Turbine Speed 
re. Design Pt 

Measured Insertion Loss (dB) 

Config. 2, 
treated nozzle, 
hardwall plug 

Config. 3, 
hardwall nozzle, 
treated plug 

Config. 4, 
treated nozzle, 
treated plug 

56% 

4.0 

0.5 

6.4 

62% 

4.6 

0.9 

6.5 

68% 

4.5 

1.3 

6.9 


Table 1: Measured insertion loss 
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III. Description of Input Parameters for FEM and CDUCT Codes 

Estimated Mode Content 

Reference [7] describes an inverse method for identifying the acoustic modal structure in 
a turbofan inlet using far-field noise measurements and high-fidelity, aeroacoustic 
radiation simulations. For the present study, that method was adapted for turbine exhaust 
noise, and used to estimate the mode content for Configuration 1 (hardwall nozzle, 
hardwall plug) at the 62% speed condition. The aeroacoustic radiation simulations were 
done using the Eversman aft radiation code. The assumed mode set included all rotor- 
stator interaction modes that could contribute to the tone at EO = 128. The mode 
amplitudes were determined by “best-fitting” the radiation calculations to the measured 
turbine-noise SPL directivity shown in Figure 6. The resultant velocity potential 
amplitudes are shown in Table 2 for the 15 largest contributing modes. This estimated 
mode content was used in the subsequent FEM and CDUCT analyses. 



Figure 6: SPL directivity pattern used to estimate mode content for Configuration 1 at 

62% design speed. 


Circumferential 
Mode Order 

Rad 

ial Mode Order 

n = 0 

n = 1 

n = 2 

n = 3 

n = 4 

m = 8 

8.82E-05 

2.39E-06 

7.50E-05 

8.57E-05 

6.40E-04 

m = -4 

7.95E-06 

2.41 E-04 

3.75E-04 

9.66E-05 

1 .92E-05 

m = 11 

2.99E-04 

2.71 E-04 

6.76E-04 

6.40E-04 

1 .88E-04 


Table 2 : Estimated velocity potential amplitudes for Configuration 1 at 62% design 

speed, EO = 128. 


Estimated Liner Impedance 

The impedance values of the hot nozzle and plug liners, Figure 2, were predicted using 
the Goodrich perforate liner impedance model [8], The liner impedance is defined as the 
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ratio of the local acoustic pressure “P” acting on the surface of the liner to the associated 
particle velocity “V” of the fluid normal to the surface (Pn & Vn): 

Zn(Normal Impedance) = Pn / Vn = (Pi+Pr) / (Vi+Vr) 

where Pi & Vi are incoming wave and Pr & Vr are reflecting wave acoustic pressure 
and particle velocity normal to the surface. Liner impedance is a function of sound 
pressure, liner geometrical dimensions, grazing flow Mach number, and boundary layer 
thickness. It is not dependent on the grazing flow direction. The general equation can be 
expressed as follows: 

Z/pc =R+jX=R 0 +R 0 f+S r Vp+R cm (V cm )+j [X m +S m Vp+X em (V cm )-cot(kh)] 

Where, 

Z/pc is a complex number representing normalized Impedance 
R is normalized acoustic resistance; X is normalized acoustic reactance 
j is V-l the unit imaginary number; p is air density and c is the speed of sound 
pc is defined as a characteristic impedance (unit: cgs- Rayl) 

R o is non-frequency dependent linear acoustic resistance 

R 0 f is frequency dependent linear acoustic resistance 
S r is the non-linear resistance slope 

is the root-mean-square particle velocity over the entire frequency range in cm/sec 
^cm is Mach Number 

Rcm(Vcm) is the grazing flow induced acoustic resistance 

is mass reactance (including ending correction) 

Sm is the non-linear mass reactance slope 
X etn (V C m) is non-linear mass reactance 
k is the wave number per cm; 

cot(kh) is the backing cavity reactance, h is cavity depth (cm) 

Pratt & Whitney provided the flow-field information (pressure, temperature, Mach 
number and displacement thickness) and internal sound-pressure spectra that were used 
as input to the Goodrich impedance model. The flow-field information was based on a 
quasi- ID flow calculation that was coupled to a boundary-layer calculation. The internal 
sound-pressure spectra were based on experimental data acquired in the turbine exhaust 
duct. 

Table 3 presents the computed impedances for the nozzle and the plug liner. 


Liner 

Location 

56% Design Speed 

62% Design Speed 

68% Design Speed 

R 

(RhoC) 

X 

(RhoC) 

R 

(RhoC) 

X 

(RhoC) 

R 

(RhoC) 

X 

(RhoC) 

Nozzle 

0.75 

-0.55 

0.85 

-0.63 

0.97 

-0.72 

Plug 

0.64 

-0.67 

0.72 

-0.75 

0.83 

-0.84 


Table 3 : Computed liner impedance at EO = 128. 
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IV. Finite Element Method (FEM) Propagation Code Predictions 

NASA and Goodrich along with University Missouri-Rolla structured a program under 
NASA contract NAS 1-00088 to study the deviations of the prediction results using 
different modal amplitudes and random phase inputs. The study results show that the 
modal phase inputs would cause significant differences in noise attenuation predictions 
[3]. G. Zlavog and W. Eversman, also reported the details of the study [1]. In this study, 
an axi-symmetric geometry model based on the approach described in References [1] and 
[3] was constructed using the modal contents shown in Table 2 as the noise source. 
These modal contents in terms of acoustic potential are related to the modal pressure 
amplitudes by: 


P = ~P o 


\dt ; 




+ U 




d£ 
dx j 


The finite element model used in the analysis is illustrated in Figure 7. The steady flow 
velocity was calculated based on the steady flow velocity potential [1] using the mean 
flow field provided by P&W at the noise source plane. The calculations were done for a 
single tone at EO = 128. The relative phase of each mode was assigned randomly, and 
1 000 sets of random phases were used for each configuration. 



Figure 7: Illustration of Finite Element Model 

Figures 8, 9 and 10 show the predicted noise attenuation levels for the treated nozzle and 
hardwall plug case, the hardwall nozzle and treated plug case, and the treated nozzle and 
treated plug case. The results are presented as histograms with each bin indicating a range 
of predicted attenuation levels. The ordinate shows the number of predicted attenuations 
that occur within each bin. As noted earlier, the total number of occurrences is 1000 for 
each configuration (i.e. 1000 sets of random-phase assignments). 

The histograms clearly show that the phase assumption significantly impacts the 
predicted attenuation values. For example, the histogram for Configuration 4 (treated 
nozzle, treated plug) shows that the predicted attenuation ranges from 3.2 dB to 13.7 dB 
depending on the mode assumption. Similar results are observed for the other two 
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configurations. Also, the histogram distributions are not Gaussian, and have positive 
skew for each configuration (i.e. they have an extended tail on the right side of the 
distribution). 


Occurences of Attenuation Based on Random Phases 
Treated Nozzle and Hardwall Plug 
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Figure 8: Predicted noise attenuation levels using random phase inputs for the treated 
nozzle and hardwall plug case, 62% design speed, EO = 128. 
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Figure 9: Predicted noise attenuation levels using random phase inputs for the hardwall 
nozzle and treated plug case, 62% design speed, EO = 128. 
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Figure 10: Predicted noise attenuation levels using random phase inputs for the treated 
nozzle and treated plug case, 62% design speed, EO = 128. 

Table 4 compares the mean predicted attenuation (i.e. the mean of the histograms) to the 
measured attenuation for each configuration. In comparing predictions to data, it is 
assumed that the predicted attenuation at EO = 128 is representative of the attenuation of 
the noise produced by all turbine stages combined. This assumption is reasonable 
because the turbine blade counts are similar for each turbine stage (i.e. between 126 and 
130). 


Configuration 

Nozzle 

Plug 

Mean Predicted 
Attenuation 
(dB) 

Measured 

Attenuation 

(dB) 

2 

Treated 

Hardwall 

4.2 

4.6 

3 

Hardwall 

Treated 

3.0 

0.9 

4 

Treated 

T reated 

7.3 

6.5 


Table 4 : LPT predicted and measured sound power attenuations at 62% design speed. 

For Configurations 2 and 4, which both have a treated nozzle, the mean predicted 
attenuations agree well with measured data (within 0.8 dB) However, for Configuration 
3, which has a hardwall nozzle and treated plug, the FEM code over predicted the mean 
attenuation by 2. 1 dB compared to the test data. Moreover, based on the histogram for 
Configuration 3 (Figure 9), the entire range of predicted attenuations exceeds the 
measured value. This result indicates that the phase assumption is not the root cause for 
the discrepancy between the predicted and measured value of attenuation. 


9 


The discrepancy between the predicted and measured attenuation for Configuration 3 
may be due to uncertainty in the plug impedance value. As described in Section III, the 
flow- fie Id information used to compute impedance was based on a quasi- ID flow 
calculation coupled to a boundary-layer calculation. Subsequent viscous CFD 
calculations have shown that the boundary-layer displacement thickness on the plug is 
much larger than that predicted from the boundary-layer calculation. Consequently, the 
actual plug resistance is probably much less than the value used in the FEM calculation. 

V. CDUCT-LaRC Predictions 

The CDUCT-LaRC [2] duct propagation and acoustic radiation tool was developed to 
study the internal acoustic propagation within and radiation from complex three- 
dimensional duct geometries. This tool was used to predict the attenuation achieved with 
the previously described acoustic treatment. The three-dimensional turbine exhaust duct 
computational grid was created using the CDUCT-LARC meshing capability. Figure 1 1 
shows the mesh between the hot nozzle and the plug. 

The mean flow used for the presented predictions was obtained from P&W viscous flow 
CFD calculations for the 68% speed condition. This solution represented a 
circumferential segment of the tailpipe geometry. The CDUCT-LaRC predictions were 
performed over the complete 3-D tailpipe geometry. Therefore, the CFD solution was 
circumferentially averaged and then interpolated onto the CDUCT-LaRC acoustic grid. 
Both the CFD and acoustic grids incorporated a fictitious plug aft of the tailplug 
termination to facilitate boundary condition specification. The transition between the 
tailplug and fictitious plug was smoothed in the CFD grid, whereas a step transition was 
present in the acoustic grid. While this had minimal effect on the mean flow 
interpolation, the geometric feature could possibly have some effect on propagation 
predictions. It should also be pointed out that the CFD solution in this case was the result 
of viscous calculations. While this is easily incorporated into the CDUCT-LaRC 
propagation calculations, it does not agree with the inviscid mean flow assumptions of 
the current CDUCT-LaRC propagation module. In this case however, the mean flow did 
not include large regions of separated flow and comparison of these predicted 
attenuations with those incorporating mean flows based on 1-D Mach-area solutions 
showed minimal changes. 

The current propagation module is based on a parabolic approximation to the convected 
Helmholtz equation, which affords predictions for complex 3-D geometries to be handled 
with relatively low computational costs. This efficiency comes at the expense of reduced 
accuracy as modes closer to cut-off are considered and the propagation angle diverges 
from the preferred angle of the parabolic approximation. Therefore, the capability to 
utilize a wide angle correction is also included in the module. The use of the wide angle 
correction can alleviate some of the error due to wide angle propagation, coming at the 
expense of increased computation time. 

Noise propagation of a combination of the first 5 radial modes for m=8, m=-4 and m=l 1 
(see Table 2) was analyzed using the CDUCT-LaRC prediction tool. In contrast to the 
FEM predictions discussed in Section IV, where random input phases were used, the 
relative phase of the various source modes was set to 0.0 and the amplitudes were set to 
those values provided in Table 2. All calculations were done for a single tone at EO=128. 
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Figure 11: CDUCT turbine exhaust duct mesh. 

Figure 12 shows the predicted Sound Pressure Levels for an untreated turbine exhaust 
duct and for the three treated configurations described in Section II. Table 5 shows the 
sound power attenuations obtained using CDUCT-LaRC along with the test results. As 
mentioned in Section IV, when comparing predictions to data, it is assumed that the 
predicted attenuation at EO = 128 is representative of the attenuation of the noise 
produced by all turbine stages combined. 

Even though the CDUCT-LaRC calculations were run with zero-phase between modes, 
the predicted attenuations agree reasonably well with the data for Configuration 2 and 4. 
Furthermore, the CDUCT-LaRC predictions with WAA are within 0.5 dB of the FEM 
random-phase predictions. These results are somewhat surprising, and may be fortuitous, 
given the phase sensitivity that was observed in the FEM calculations. The wide-angle 
approximation (WAA) increased the predicted attenuation, which improved agreement 
with test data for Configuration 2, but worsened the agreement for Configurations 3 and 
4. Like the FEM code, CDUCT-LaRC over predicts the attenuation for Configuration 3. 
This discrepancy may be due to uncertainty in the plug impedance (see discussion in 
Section IV). 

A random-phase CDUCT-LaRC study was also conducted similar to the FEM study 
described in Section IV. The flow field, mode inputs and impedances were the same as 
previously described, but 1000 cases were run with a random phase assigned to each 
mode. The wide-angle approximation was not used during the CDUCT-LaRC random- 
phase study. 

Figure 13 shows a histogram of the CDUCT-LaRC attenuation results for Configuration 
4 (treated nozzle, treated plug). The predicted mean (6.1 dB) compares well with the 
measured attenuation (6.5 dB), and is close to the zero-phase result (6.6 dB). The 
CDUCT-LaRC histogram also agrees well with the FEM code, showing a broad range of 
attenuation with positive skew. 
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Hardwall Nozzle 



Hardwall Plug 




Hardwall Plug 
Treated Nozzle 



Treated Plug 
Hardwall Nozzle 



Figure 12: SPL predictions at EO = 128 for four configurations, 62% design speed. 
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Config 

Nozzle 

Plug 

Percent 

Design 

Speed 

Predicted 

Attenuation 

(dB) 

No WAA 

Predicted 

Attenuation 

(dB) 

With WAA 

Measured 

Attenuation 

(dB) 

2 

Treated 

Hardwall 

56 

2.9 

4.3 

4.0 

62 

2.9 

4.3 

4.6 

68 

2.9 

4.3 

4.5 

3 

Hardwall 

Treated 

56 

2.8 

3.3 

0.5 

62 

2.9 

3.5 

0.9 

68 

3.1 

3.6 

1.3 

4 

Treated 

Treated 

56 

6.4 

7.6 

6.4 

62 

6.6 

7.8 

6.5 

68 

6.9 

7.9 

6.9 


Table 5 : CDUCT-LaRC predicted and measured sound power attenuations. Predictions 
with and without Wide Angle Approximation (WAA). 


TREATED NOZZLE AND TREATED PLUG 



Figure 13: CDUCT-LaRC predicted attenuation for treated plug and nozzle, random 

phase, 62% design speed, EO = 128 

VI. Summary 

NASA, Pratt & Whitney and Goodrich have validated an approach for predicting turbine 
noise attenuation on a full-scale engine. The approach consists of three steps: (1) 
estimating the turbine-noise mode content, (2) computing the liner impedance based on 
the computed flow-field, and (3) computing in-duct propagation using high-fidelity codes 
such as the Eversman FEM code and CDUCT-LaRC. The predicted noise attenuation 
agreed to within 1 dB of the measured attenuation, except for one case. In that case, the 
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discrepancy between the predicted and measured attenuation could be the result of 
uncertainty in the plug impedance. Further studies are needed to confirm this hypothesis. 

In the noise propagation analysis, the modal phase assumption significantly impacts the 
predicted attenuation values. However, when averaged over a large number of random- 
phase assumptions, the predicted attenuations are in good agreement with the data. 
Therefore, the recommended approach is to randomly assign phase to each mode, run a 
multitude of random-phase assumptions, and use the mean value to represent the 
predicted attenuation. 

With the wide-angle approximation, the CDUCT-LaRC results were within 0.5 dB of the 
Eversman FEM results. This result supports the use of CDUCT-LaRC for the predicting 
the noise attenuation of more complex 3D geometries (e.g. fan ducts with pylon and 
bifurcations) that cannot be handled with the 2D Eversman FEM code. Further studies are 
needed to validate the approach for more complex geometries on a full-scale engine. 
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